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ABSTRACT 

We present multi-wavelength, high spatial resolution imaging of the IRS 7 
region in the R Corona Australis molecular cloud. Our observations include 
1.1 mm continuum and HCO+ J = 3 ^ 2 images from the SMA, ^^CO J = 3 ^ 2 
outflow maps from the DesertStar heterodyne array receiver on the HHT, 450 /im 
and 850 /im continuum images from SCUBA, and archival Spitzer IRAC and 
MIPS 24 nm images. The accurate astrometry of the IRAC images allow us to 
identify IRS 7 with the cm source VLA lOW (IRS 7A) and the X-ray source Xw- 
The SMA 1.1 mm image reveals two compact continuum sources which are also 
distinguishable at 450 /im. SMA 1 coincides with X-ray source CXOU J190156.4- 
365728 and VLA cm source lOE (IRS 7B) and is seen in the IRAC and MIPS 
images. SMA 2 has no infrared counterpart but coincides with cm source VLA 9. 
Spectral energy distributions constructed from SMA, SCUBA and Spitzer data 
yield bolometric temperatures of 83 K for SMA 1 and <70 K for SMA 2. These 
temperatures along with the submillimeter to total luminosity ratios indicate that 
SMA 2 is a Class protostar, while SMA 1 is a Class 0/Class I transitional object 
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(L=17±6 L0). The ^^CO J = 3 ^ 2 outflow map shows one major and possibly 
several smaller outflows centered on the IRS 7 region, with masses and energetics 
consistent with previous work. We identify the Class source SMA2/VLA9 as 
the main driver of this outflow. The complex and clumpy spatial and velocity 
distribution of the HCO^ J = 3 ^ 2 emission is not consistent with either bulk 
rotation, or any known molecular outflow activity. 

Subject headings: stars: formation, ISM: individual - R CrA - Clouds - ISM: jets 
and outflows - ISM: kinematics and dynamics - ISM: molecules - radio lines: ISM 



1. Introduction 



In the past 15 ye ars, the search fo r Class protostellar objects has yielded only ~100 
unambiguous sources (jFroebrich I l2005l ) . Class sources, which display a high ratio of sub- 
millimeter to bolometric luminosity, ofte n remain undetect able in the near-IR, and show 
bolometric temperatures of 80 K or less (IMyers et al.lll998l ). An assortment of recent ob- 
servations concentrating on the core of the R Corona Australis molecular cloud (R CrA) 
has revealed sources that could be Class objects. R CrA is an ideal laboratory for the 
study of low mass star fo rmation because at a distance of ~170 pc, it is one of the nearest 
examples of this activity (IKnude fc H0g||l998l ). As a result, observations have high physical 
resolution (170 AU per arcsecond) when compared to more distant star forming regions. 
While images with the Chandra X-ray Observatory and the VLA have been able to resolve 
individual sources in the R CrA region, work in the (sub)millimeter band has been limited 
to single dish telescopes, with > 20" resolution. Nevertheless, the proximity of R CrA has 
allowed even moderate resoluti on studies to reveal (sub)millimeter sources in the region th at 
have no near-IR counterparts (jChini et al.ll2003l : iGroppi et al.ll2004al : ISchoier et al.ll2006l ). 



Early single-dish observations of this region wit h a 2.3 ^ beam identifled a molecular 
cloud core containing the emission line star R CrA (jhoren 1979 ) . Later observati ons of 
C^^O with a 45" beam re solved this core i nto na ultiple components (iHarju et al.lll993l ). Ob- 
serving with a 23" beam, iHenning et al. I (jl994j ) attributed the bulk of th e long wavelength 



dust continuum emission to the deeply embedded near-IR source IRS 7 (ITaylor fc Storey 



19841). A raolecular outflow approximately centered on IR S 7 has been known for many years 
Levreault Ill988l : iHarju et al.lll993l : iGroppi et al.ll2004al. and others), and is an indicator of 



protostellar activity. More recent work by lNutter et al. 



(120051 ) with the JCMT at 850 /im and 



450 ixm has begun to disentangle the submillimeter emission into numerous sources in the 
vicinity of IRS 7. It is now clear that IRS 7 is unlikely to be the powering sourc e of this long 



wavelength dust emission. The submillimeter sources SMMIB and SMMIC (iNutter et al. 
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20051 ) are associated with two centimeter continuum point sources (IFeigelson et al 



These centimeter sources are also detected in the X-ray by iHamaguchi et aLl (120051 ). with 



19981). 



properties suggestive of Class or very young Class I objects. Clearly, millimeter or sub- 
millimeter interferometric observations are necessary to unambiguously identify the young, 
embedded protostellar sources in the region. 

We have observed the central region of the R CrA molecular cloud with the Submillime- 
ter Array on Mauna Kea, Hawaii in both continuum (1.1 mm) and spectroscopic 

(HCO^ J = 3 — > 2) mode, with the goal of spatially resolving and characterizing the Class 
and Class I protostellar candidates in the region. In this paper, we combine the high spatial 
resolution SMA images with archival Spitzer IRAC and MIPS photometry, JCMT 850 /im 
and 450 fim maps and ^^CO J = 3 ^ 2 outflow maps of moderate spatial resolution in 
order to more clearly reveal two young protostellar objects in the vicinity of IRS 7. Our 
observations range in linear resolution from 0.02 pc to 0.001 pc. 



2. Observations 

2.1. SMA 1.1 mm Continuum and HCO+ J = 3 ^ 2 Observations 

R CrA IRS 7 (asooo.o = 01"* 55.^^34, (52ooo.o = -36° 57' 21^7) was observed with the 
SMA on two occasions. The first track was taken on May 23, 2004 with a somewhat extended 
array configuration, with pads 1, 7, 8, 11, and 23 used. The zenith opacity at 225 GHz as 
measured with the tipping radiometer at the Caltech Submillimeter Observatory (CSO) 
ranged from 0.14 to 0.20 during the observations. The tuning was chosen to place the HCO"*" 
J = 3 ^ 2 line (267.557619 GHz) in correlator chunk 3 of the lower sideband with f lsr = +6 
km s~^ (corresponding to a local oscillator setting of 271.77 GHz). The spectral resolution 
was 256 channels per correlator chunk (0.312 MHz, or 0.35 km s~^ per channel). Uranus and 
Neptune were observed for bandpass and flux calibration, while the quasar J1924-292 (5.5 
Jy) was used as the primary gain calibrator and pointing source. The UV data span 8 to 100 
kX, with the greatest coverage between ~20 and 60 kX. The primary beamsize was 46". The 
second track of data was taken on October 4, 2004 in the compact-north configuration, with 
antennas on pads 1, 5, 7, 8, 9, 10, 12 and 16. Seven of the eight antennas provided usable 
data for this track. These data were taken in a hybrid resolution mode with 512 channels 
per chunk (0.156 MHz, 0.18 km s~^ ) for the HCO"'" J = 3 — >^ 2 correlator chunk, and 128 



^The Submillimeter Array (SMA) is a collaborative project between the Smithsonian Astrophysical Ob- 
servatory and the Academia Sinica Institute of Astronomy & Astrophysics of Taiwan. 
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channels (0.625 MHz, 0.70 km s~^) for all other chunks. In an attempt to detect the source 
SMM lA, a second source position was observed one half of a primary beam south of IRS 7. 
The calibration sources and observing strategies were otherwise identical to the first track. 

These data were calibrated and reduced with the MIRIAD sof tware package with ex- 



tensions for the SMA using the procedures outlined in IZhao I (120061 ). The continuum data 
were constructed using all the line-free "chunks" of the correlator to produce images for 
the north and south pointings. The south pointing was not added into the final image, 
because no source was detected. The continuum data were inverted to the image domain 
using uniform weighting, then cleaned using a gain of 0.08 with 2500 iterations. An rms 
noise of 9 mJy/beam was achieved. The synthesized beam is 4''2 x 1'.'7 (700 x 300 AU) with 
a position angle (p. a.) of 27°. We estimate the absolute astrometric uncertainty to be 0''3 
and the flux scale uncertainty to be 20%. The spectrometer "chunk" with the HCO^ J = 
3 — 2 line was extracted from the data of the two tracks and independently processed. The 
data from the second track were binned to match the spectral resolution of the first track 
prior to processing. The data were then calibrated and reduced in a similar fashion to the 
continuum data. Inversion to the image domain was performed with natural weighting for 
maximum sensitivity. Cleaning was performed with a gain of 0.08 and 2500 iterations. As 
with the continuum, no line emission was detected in the southern pointing. 



2.2. Heinrich Hertz Telescope CO J = 3 ^ 2 Observations 

The ^^CO J = 3 ^ 2 data for this work was taken with the 10 m Heinrich Hertz Telescope 
at Mt. Graham, AZ d uring April, 2006 with the DesertStar 7 beam heterodyne array receiver 



(iGroppi et al.ll2004bl ). 5' x 5' maps centered on IRS 7 were taken in the on-the-fly (OTF) 
observing mode, using the facility 1 MHz resolution filter bank spectrometers. Eight OTF 
maps were made over the course of the run, and later coadded. The spectral resolution is 
0.87 km s~^, with a spatial resolution of 23". Hot-sky calibrations were done once per row, 
with cold load calibrations done once per map. Data were gridded and combined based 
on RMS weighting using the CLASS software package. The intensities of each pixel were 
scaled using their measured main beam efficiency before combining. Further analysis was 
performed using the IDL package. 
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2.3. James Clerk Maxwell Telescope 850 /im and 450 fim Observations 

The 850 fim and 450 /im continuum observat ions were obtained witli bolometer array 
SCUBA on JCMtI, Mauna Kea, Hawaii. SCUBA jHoUand et al.lll999h has 37 bolometers in 



the long and 91 in the short wavelength array separated by approximately two beam widths 
in a hexagonal pattern. The field of view of both arrays is ~2'3. Both arrays can be used 
simultaneously by means of a dichroic beamsplitter. 

The SCUBA ob servations reported in this paper were all obtained in jiggle-map mode 



( iHoUand et al.lll999l ) with chop throws of 120" - 150" in azimuth. Because the R CrA cloud 
complex is rather extended, it was impossible to completely chop off the array, but any 
strong emission in the off source positions was carefully blanked out in the data reduction 
stage and should not affect the morphology or photometry of the point-like sources in the 
field. It does, however, affect the integrated flux density of the extended cloud emission, 
which is irrelevant for our analysis, since we are only interested in the embedded compact 
sources in the cloud. 

The data used in this paper come from three 5- integration maps obtained during SCUBA 
commissioning on March 30 and April 5, 1997. These maps were obtained in excellent submil- 
lim eter conditions w i th CS Q r225GHz ~ 0.045 - 0.05. Analysis of this data set was discussed 
by van den Ancker ( 1999 ). We additionally searched the JCMT archive at the CADC§| 



and found one more 6-integration map3 obtained in good submillimeter conditions (CSO 
T225GHZ ~ 0.066) on April 11, 2000, which we added to our data set. Pointing corrections 
were derived from observations of the blazar 1929-293 and flux calibration from observations 
of Uranus. We estimate the calibration accuracy to be better than 10% at 850 fim and 
better than 20% and 450 yum. The measured Half Power Beam Width (HPBW) for the 1997 
and 2000 observations are remarkably similar. In 1997 we used a 120" chop throw, but an 
experimental chop waveform, which broadened the beam in the chop direction, while the 
observation from April 2000 used a 150" chop throw, which also broadened the beam by 
about the same amount. The HPBW is therefore ~ 15"6 x 13"5 and 9"1 x 7"8 for 850 fim 
and 450 /xm, respectively with the major axis aligned with the chop direction. 



^The JCMT is operated by the Joint Astronomy Centre, on behalf of the UK Particle Physics and As- 
tronomy Research Council, the Netherlands Organization for Scientific Research, and the Canadian National 
Research Council. 

•^Guest User, Canadian Astronomy Data Center, which is operated by the Dominion Astrophysical Ob- 
servatory for the National Research Council of Canada's Herzberg Institute of Astrophysics 

''This map is from the data set published by Nutter et al. ( 2005[ l 
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The data were reduced in a standard way using SURF (jjenness and Lightfoot I Il999 



Sandell. Jessop and JennessI l200ll ) and STARLINK imaging software, i.e., we flat fielded, 
extinction corrected, sky subtracted, despiked, and calibrated the images in Jy beam~^. 
We then pointing corrected each data set for any drift in pointing between each pointing 
observation and added the data together to determine the most likely submillimeter position 
at 850 fj,m. Once we had derived a basic 850 fim astrometric image, we made Gaussian fits 
of point-like sources for each data set and derived small additional RA and Dec corrections 
to each scan (shift and add) to sharpen the final image to this position. Manual Gaussian 
fitting for these closely spaced, embedded sources is more reliable than automatic clump 
finding routines (e.g. clumpfind), which often to not find clumps embedded in an extended 
cloud, especially if those clumps are close to each other. Since there is a small mis-alignment 
between the 850 and 450 yum arrays, we first corrected the 450 /xm images for any pointing 
drifts and then did shift and add using an isolated compact source in the 850 /im image as 
an astrometric reference source for the 450 /im data. We estimate the astrometric accuracy 
in our SCUBA images to be < 2". 

The final coadd was done by noise- weighting the data in order to minimize the noise in 
the final images. The images were coadded by weighting with the inverse square of the noise 
in each set and accounting for the difference in integration time. The first image is taken as 
a reference and then the weighting factors are computed for each subsequent image: 



W2 



t2/h Y 



This method minimizes the noise in the final coadded image. The rms of the 450 /xm im- 
age is ~ 0.3 Jy beam~^ and ~ 12 mJy beam~^ f or the 850 /xni image. All our ma ps were 
converted to FITS-files and exported to MIRIAD (jSault. Teuben and Wright Ill995l ) for fur- 
ther analysis. In order to correct for the error lobe contribution, especially at 450 fim, we 
have deconvolved all the maps using CLEAN and an elliptical model beam derived from our 
Uranus observations. Our model beam is composed of three symmetric Gaussians: a main 
beam, a near-error lobe, and an extended, low amplitude, far-error lobe. We rotated this 
model beam to the approximate rotation angle of each individual map and took the average 
of the rotated beams to provide a more accurate model beam for our coadded maps. 
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3. Analysis 



3.1. SMA continuum results 



Figure [T] shows the SMA 1.1 mm continuum image. We detect two sources of approx- 
imately equal flux density. The positions and flux densities derived from two-dimensional 
Gaussian flts to the sources are given in Tabled! SMAl, the eastern source, is marginally 
resolved, with an estimated source size of l''l x 0''5 at p. a. = 58°. The western source, 
SMA 2, is unresolved. SMAl coincides to within 0.^1 6 of a centi meter source designat ed by 
various names: V LA lOE (IRS 7B) by 
source 14/IRS7E JPorbrich et al.ll2006 



Brownl (Il987l ). source 8 (IFeigelson et al.lll998l ). and 



2007). SMA 1 also coincides within 0^'17 of the hard 



X-ray source CXOU J190156.4-365728 (IHamaguchi et al.ll2005f ). SMAl is detected in all 
IRAC bands and at 24 fim with MIPS (see Section [3.31 and Table [3D • In contrast, SMA 2 
is not seen in either the X-ray, near-IR or mid-IR bands, but it s position is co nsistent with 
a centimeter source d esignated by vari ous names: sou r ce 9 by iBrownl (119871 ) , source 6 by 
Feigelson et aD Jl998[ ). and source 11 bv lForbrich et al.l J2006h. The SMA position of SMA 2 
is within 0''5 of the S-band p osition reporte d by Brown (jl987 ) which agrees to O'.'l of the 
X-band contour map peak of iForbrich et al.l (120061 ). Although those authors do not quote 



position uncertainties, IFeigelson et al.l (Il998l ) notes that the VLA position accuracy for their 
data on this low declination (i.e. very low elevation) source is ±1". Even if the VLA position 
is accurate to ±0.2", this uncertainty combined with the SMA position uncertainty allow 
for the possibility that the 1.1 mm emission from SMA 2 is coin cident with the centimete r 
source. SMA 2 also is within 0('93 of the 7 mm VLA source CT2 Jchoi fc Tatematsu 112004 ): 
the implications of this possible association are disc ussed later. SMA 1 and SMA 2 coincide 
with the two sub-mm sources SMM IB and SMM IC lNutter et al.l (120051 ) . which is seen more 
clearly in our SCUBA data (Section [S^D- We do not detect IRS 7A (VLA lOW, Feigelson 7, 
Forbrich 12/IRS7W), nor do we detect any sources in the "prestellar" core SMM lA, which is 
the brightest extended source in the submillimeter continuum, located ~20" south of IRS 7. 



3.2. SCUBA 450 /xm and 850 /xm results 



Our 850 Atm and 450 fim S CUBA images (Figure [2]) agree qualitatively with the SCUBA 
results of iNutter et al.l (120051 ) . although we have better signal-to- noise and image quality 
due to a larger data set and imp r oved data processing methods. We detect all the sub- 
mm sources seen by iNutter et al.l (120051 ) and additionally we identify a fainter core south 
of SMM lA. Although SMM lA is not resolved at 850 /im, the 450 fim image suggests that 
it consists of two cores, rather than a single elongated elliptical core. This was also seen 
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by Ivan den Ancker I (119991 ) in a high spatial resolution 800 /im-map and Ivan den Ancker 
therefore listed SMM lA as two separate sources. We have used the MIRIAD task IMFIT 
to fit elliptical Gaussians with a two-component Gaussian fit; one compact component for 
the core and an extended Gaussian for the surrounding cloud. Even though the surrounding 
cloud may not be well described by a single Gaussian, this approach works reasonably well if 
one limits the region over which the fitting is done, so that the extended cloud emission can 
be approximated by a single Gaussian. For very extended cores like SMM 1 A and SMM 4 we 
approximated the surrounding cloud emission with a baseline plane rather than a Gaussian. 
The uncertainty in the estimation of the background, however, adds about a 10 - 15% error 
to the integrated flux densities. The integrated flux densities and positions for the sub- 
millimeter continuum sources ar e given in Table El These flux densities differ significantly 
from the flux densities quoted by lNutter et al.l (120051 ) . because their results include emission 
from the surrounding cloud. 

Only two of the known near-IR sources (IRS 1 and WBM 55) in the R CrA cloud core 
coincide with compact submillimeter SCUBA sources, but in our high spatial resolution 450 
fim image we also resolve the two 1.1mm SMA continuum sources, SMA 1, and SMA 2 as two 
separate compact dust emission peaks. IRS 5 and the new mid-IR source IRAC 9 (Table 3) 
lie close to the peak of the extended submillimeter core SMM 4, but we see no detectable dust 
emission excess from the stars themselves, even though our ^^CO J = 3 — > 2 observations 
(Section 13.41) indicate that there may be a molecular outflow powered by one of these stars. 
Neither do we detect dust emission from the two y oung HAEBE stars R CrA and T CrA. 
R CrA is the most luminous source in the region (IWilking. Taylor &: Storey Ill986l ). It is 
very strong in the mid-IR and illuminates a large reflection nebula. It is possible that the 
star may still be associated with a compact dust disk, but since it is situated near the cloud 
ridge encompassing SMA 1 and 2 (Figure [2]), it is difficult to detect it, especially if the disk 
is small. From our 450 /im image we estimate that the emissi on from such a disk would be 
less than 0.7 Jy beam~^. Using the dust grain formalism of iHildebrand I (119831 ). this flux 
density limit corresponds to a mass limit of < 1.2 x 10^^ M0 (assuming T = 30K and = 1, 



i.e. the average disk SED parameters found in the survey of I Andrews fc Williams! (120051 )) 



IRS 1, also known as HH 100-IR, illumi nates a prominent cigar-shaped reflection nebula 
(IHartigan fc Graham 1119871 : iGraham Ill993l ) with an HH object (HH97) near the tip of the 
nebula, indicating that one sees the cavity walls of an outflow driven by IRS 1. WBM 55 is 
near the HH objects HH104 A & B, but these are probably excited by SMA 2 (see Section 
3.4p or some other young star in the the vicinity of SMA 2. Nevertheless it is likely that 
WBM 55 also drives an outflow. However, we do not see any dust emission from IRS 7, 
which has been suggested as a poss i ble candidate for d riving the large scale CO outflow 
Jlevreault IligSsUcroppi et al.ll2004al : Iwilking et al.lll997h . IRS 7 is deeply embedded {Ay > 
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35"^) in the cloud core and appears as a nebulous knot in K band. IWilking et al.l (119971 ) 
suggested that the 2.2 /im emission may originate from the nebulosity and that the star 
itself is invisible at K band. This sug g estion gets further support f r om th e high contrast 
images presented by IChen &: GrahamI (119931 ) and iHamaguchi et al.l (120051 ) which show a 
larger reflection nebulosity to the northeast seemingly forming a bipolar nebula. In this 
scenario one would expect the mid-IR (i.e. IRAC emission) and VLA radio source to pinpoint 



the lo cation of the star and place it northwest of the near-IR position. IHamaguchi et al. 



(120051 ) derive an accurate 2.2 /xm position, which is slightly north of the IRAC/VLA position, 
suggesting that the 2.2 /xm emission has the same energy origin as the mid-IR emission. To 
illustrate the complexity of this region we display a sub-image of the 450 /im image covering 
only SMAl, IRS 7 and SMA2 in Figure [3l Here w e overlay the IRAC positions (Table 
3), the cm VLA po s itions from iFeigelson et al.l (119981 ). and the 7 mm VLA positions from 



Choi &: Tatematsu I (120041 ). This figure shows clearly that SMA2 and IRS 7 are different 
objects. 



3.3. Spitzer IRAC and MIPS archival data 



We retrieved Spitzer GT O observations for the R CrA region from the Spitzer data 
archive (iMegeath et al.l 120041 ). These GTO observations included IRAC images at 3.6 /im, 
4.5 /xm, 5.8 /im, and 8.0 /xm, as well as MIPS observations at 24 /xn]^. Figure S] shows a 
composite IRAC 4.5 /im image derived from the "long" and "short" exposure time images 
for increased dynamic range. The spatial resolution of IRAC permitted the use of simple 
aperture photometry using the Starlink Gaia package to extract positions and flux densities 
for all 8 /im sources within a 2' radius of IRS 7 (Table [3]). The photometric uncertainty is es- 
timated at ~10%, dominated by the uncertainty of the background level in the artifact-filled 
field. The astrometric accuracy of the IRAC images (< 0''3) con firms that th e infrared object 



IRS 7 is associated with the VLA so urce IRS 7A (VLA lOW) (ISrownl 119871 : IFeigelson et al. 



1998) and the hard X-ray source X^i/ ( Hamaguchi et al. 20051 ). For the first time we also de- 
tect an infrared counterpart to the radio source IRS 7B (VLA lOE), which we have identified 
as SMA 1 in our 1 . 1 mm continuum imaging. We have also discovered a new, relatively faint 
mid-IR source, IRAC 9, ~ 9" northeast of IRS 5. 

The pipeline reduced "Post-BCD" MIPS images were used to obtain aperture and PSF 
fitted photometry for any 24 /xm sources in the vicinity of SMA 1 and SMA 2. SMA 1 was 



'''The spatial resolution of the 70 and 160 /im MIPS data were too coarse to reveal any source information 
in the R CrA cloud core. 
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clearly detected at 24 /iin but it was saturated and very near the sources IRS 7 and R CrA, 
which were also badly saturated. We estimated flux densities by fltting the unsaturated 
wings of the PSFs to create a "restored" unsaturated PSF for flux estimation. The process 
was iterative. We flrst flt R CrA (200±20 Jy), and subtracted the model PSF from the image, 
then flt IRS 7 (65±15 Jy), and flnally flt SMAl (20±5 Jy). Th e flux density obtained by 



this p rocedure agrees quite well with ground based photometry. IWilking. Taylor fc Storey 



(119861) quote 218 ± 6 Jy for R CrA and 37 Jy for IRS 7 at Q (20 /xm). After the subtraction 
of the three model PSFs, there was no indication that SMA2 was detected. The fltting 
and subtraction of the highly saturated PSFs likely masked any faint emission from SMA2. 
We use the flux density uncertainty estimate of SMA 1 to place an upper limit on the flux 
density of SMA 2 < 5 Jy. 



3.4. CO outflows 



We have constructed outflow maps from ^^CO J = 3 ^ 2 HHT data which provide 
moderate (23") spatial resolution. Figure [5] shows the CO outflow red and blue contours. 



The ve locity intervals used to produce the red and blue contours are similar to iGroppi et al. 
( l2004al ). at -10 to +3 km s^^ for the blue wing, and 10 to 21 km s~^ for the red wing. 



These data provide about three times the spatial resolution of iGroppi et al.l (l2004al ). and 
conflrm that the center of the outflow is in the vicinity of IRS 7. Even with this spatial 
resolution the CO map gives the appearance that the high velocity emission is dominated 
by a single outflow centered near IRS 7 and SMA 2. A closer inspection, however, shows 
features which are hard to explain with a single outflow. This is not surprising, since we 
know that there is a young cluster embedded in the molecular cloud core, several of which 
appear to be surrounded by a ctive accretion d isks. We plot the optica lly-detected emission 
line objects (HH objects) from lGraham I Jl993h and lWang et al. I (120041 ) in flgureO In flgure 
|5]we draw a line connecting the peaks of the red and blue outflow lobes. This line (at p. a. 
105°) passes almost exactly through SMA 2, and also is very close to the HH104 objects. 



Wang et al. I (120041 ) beheve that HH104 C and D are associated with IRS6 and a nearby 
optically detected jet, and not associated with HH104 A and B. The spatial association of 
IIII104 C and D with the SMA 2 outflow axis could therefore be coincidental. SMA 2 is also 
located at the approximate EW center of the CO outflow lobes. We believe this is sufficient 
evidence to attribute the major outflow activity in the region to SMA 2. Further evidence 
for this concl u sion comes from the fact that in all 13 of the molecular lines measured by 



Schoier et al.l (120061 ) . the FWHM linewidths are broader when the APEX telecope beam 



encompasses SMA 2, compared to SMAl. 
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We also find a second, smaller CO outflow, centered on the infrared source HHIOO-IR 
(IStrom et al.lll974j ). with the outflow axis parallel to the line connecting the HH objects 
HH98 (slightly red-shifted) and the two blue-shifted HH- obiects HHIOO and HH97. This 
outflow could also be associated with HH99A and HH99B (ICohen et al.lll984l ). A measured 
proper motion for HH99 reported by ICaratti o Garatti et all J2006h is at a similar position 
angle to the line connecting the outflow lobes shown in [51 The blue-shifted outflow lobe is 
less prominent than the optical/near-IR cavity, presumably because the blue outflow emerges 
out of the dense cloud core into the low density halo, and therefore encounters much less gas 
to interact with. Other small outflows could exist, but do not have the reinforcing evidence 
of the presence of optically-detected emission line objects. SMAl, for example, could be 
driving an outflow oriented NNW, which could explain the blue-shifted ridge northwest of 
SMA2, although there is also faint blue-shifted emission to the south of SMAl. Another 
possibility for the blue-shifted ridge NW of SMA 2 is a roughly E-W outflow from IRS 6. 
Alternatively, the red-shifted emission in the vicinity of IRS 6 could be part of a NE-oriented 
red-shifted outflow lobe originating from or near IRS 5 with the blue outflow lobe too faint 
to be detected in the low-density cloud SW of IRS 5. One will need much higher spatial 
resolution and sensitivity to disentangle the outflows from the young stars in the R CrA 
core. Therefore, in the following analysis, we compute mass estimates from the integrated 
red-shifted and blue-shifted emission as if it originates from a single outflow. We assume that 
the excitation temperature is equal to the peak CO line temperature at each map position 
(i.e. the line cores are optically thick), and the line wings are opticall y thin. We c ompu te 
the total mass and energetics for the entire outflow complex. Following ISnell et al.l (119841 ). 



N, 



CO 



1.94 X 10 kTe^ 
Aui hB^Qu 



' Tn{CO)dVexp(-^ 

^ J- P.T. 



era 



(2) 



We calculate a red lobe mass of 0.04 M©, a nd a blue lobe ma ss of 0.03 M©, assuming a 
CO abundance relative to H2 of 1.5 x 10~^ ( Hogerheiide I 2003). T his is considerably less 



than the masses calculated for the outflow in iGroppi et al.l (j2004al l. but these data cover 
only a quarter of the area on the sky. In addition, the assumption that the line wings are 
op tically thin m a ke thes e estimates lower limits for the mass. As shown for ^^CO J = 1 — 



m iGroppi et al.l (j2004al ). this assumption leads to an underestimation of the mass in the 
outflow by a factor of several, depending on the level of self absorption in the line. Also, the 
filling factor of the warmer gas traced by ^^CO J = 3 ^ 2 could be less than in the case of 
^^CO J = 1 ^ 0. With these caveats, the mass and associated m omentum, kiii e tic ene rgy 
and mechanical luminosity are consistent with the results found in IGroppi et al.l (l2004al ). 
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3.5. Interferometric HCO+ J = 3 ^ 2 Observations 

The interferometric HCO"'" J = 3 ^ 2 observations reveal a very clumpy distribution 
of emitting gas with a complex velocity field. Clearly, the lack of short spacing information 
severely limits the contents and fidelity of the image. Lines are typically about 2 km s~^ 
wide, and can exhibit complex line shapes. A simplified contour plot with "red" and "blue" 
contours is shown in Figure [6l This plot shows that the clumps of HCO"*" J = 3 ^ 2 emission 
seem to follow a roughly spherical distribution around SMA 2, and show no systematic spatial 
distribution of velocity. In addition, the HCO"*" J = 3 ^ 2 clumps are not spatially coincident 
with either of the continuum sources SMA 1 or SMA 2. The velocities of the clumps are also 
not consistent with the velocity distribution of the molecular oufiow seen in ^^CO J = 
3^2. The observed clumpy, chaotic velocity field also does not support the hypothesis of a 
circumstellar or circumbinary disk , as was suggested b ased on the kinematics of single-dish 



observations of HCO+ J = 4 ^ 3 fIGroppi et al.ll2004af ) 



This clumpy structure might be the result of the interaction of many oufiows from the 
large number of protostars in the region. The velocity field of the clumps does not seem to 
correlate with the large scale outfiows seen in CO J = 3 ^ 2, but the complexity of the 
region is high, and the interferometric image obviously suffers from missing flux. Another 
explanation for the observed morphology is that the clumps are knots in a spherical shell 
of material swept-up by the molecular outflows and winds produced by the young cluster. 
The critical density of the HCO+ J = 3 ^ 2 transition in which the clumps are observed 
suggests the density is of the order ~ 10^ cm^'^. In addition, the spatial flltering provided 
by the interferometer also favors the detection of clumpy structure. If this is the case, we 
can calculate an upper limit on the mass of the possible gaseous shell by assuming the gas 
is distributed in a uniform sphere, and that the column density is that of the brightest 
clump. Assuming the lines are optically thin, we can calculate the column density in the 
same manner as the previous section. We can then estimate the mass by integrating over a 
circle of diameter 45 " (8.5 x 10^^ cm at a distance of 170 pc). We convert the flux density 
to brightness temperature using the following formula: 

TMK) ^ , "■"^^-''^g..... (3) 



, I 3.92x10-8 u{GHz)3 gife ^ 



where 9i and 62 are the FWHM dimensions of the source region (in this case, the beam size). 
This results in a HCO"*" J = 3 ^ 2 brightness temperature of 17.7 K, and an i ntegrated in- 



tensi ty of 24.3 K km s ^ . Assuming an abundance of HCO"*" to H2 of 1.2 x 10 * (jHogerheijde 
2OO3I ) ■ we calculate a shell mass of 2.1 x 10'^^ g, or 0.01 M©. With an expansion velocity of 
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~ 1.5 km s~^, the momentum and kinetic energy of the shell is small compared to that of the 
large scale molecular outflow. Thus, the driving source(s) of the outfiow(s) are easily capable 
of driving the expansion of such a shell, although the mechanism is unclear. Alternatively, 
the shell could be driven by a wind from a more evolved pre-main-sequence star like IRS 7. 

It is important to note that the ring like morphology we see in the interferometric 
HCO"^ map could be the result of a lack of zero spacing and short spacing data, rather than 
a reflection of the true gas morphology. Future observations with combined interferometric 
and single dish data will be necessary to generate an unambiguous determination of the gas 
morphology around SMA 2. 



4. Discussion 

4.1. FIR Spectral Energy Distributions 

The SMA 1.1mm data and SCUBA 850 fim and 450 /zm data are combined with the 
MIPS and IRAC archival data to produce SEDs capable of constraining the peak of a grey- 
body, thereby estimating the source temperature. The source SEDs and greybody fits are 
shown in Figures 7| and El Neither sou rce is detected in the near-infrared. The upper limits 
plotted are from IWilking et al.l (119971 ) who quote 3a upper limits of (17.5, 17.0 16.5) for 



the (J, H, K) bands. We fit SMA 1 with a two-temperature composite greybody. The 
uncertainties in the fitted and derived parameters have been estimated by perturbing the 
flux densities by la (in such a way as to maximize the change in fitted temperature) and 
repeating the fit. We find temperatures of 55 ± 3 K and 344 ± 6 K for the two components, 
with a low frequency dust grain emissivity index (3 = 1.3 ± 0.2. The warmer greybody 
was assumed to have (3 = 1.0. The integrated bolometric luminosity based on the SED for 
SMAl is 17 ± 6 Lq. The total mass is (3.5 ± 1.4) x 10^^ Mq, assuming a gas to dust 
mass ratio of 100. The ratio of submillimeter luminosity to total luminosity in SMA 1 is 
(1.4 ± 0.4) X 10~^, suggestive of a class O/I source. The near infrared flux density upper 
limits are a fac tor of 2.5 x 10^ less than the SED peak 260 Jy at 62yum), suggestive of a 



Class source (Froebrich 



2005*). We also calculate the bolometric temperature of the source 
following iMyers fc Ladd I (Il993i) . The bolometric temperature T^oi of a SED F^, is defined as 



the temperature of a blackbody having the same mean frequency u: 



Tbo; = [C(4)/4C(5)]/iz>A = 1.25 x 10 



"z/ K Hz-\ 



(4) 



where ({n) is the Riemann zeta function of argument n, h is Planck's constant, k is Boltz- 
mann's constant, and the mean frequency u is the ratio of the first and zeroeth frequency 
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moments of the spectrum: 



Jo 



(5) 



We numerically integrate the fitted SED to obtain a bolometric temperature of 83 K for 
SMAl. 

For SMA2, only an upper limit exists for the 24 yum flux density, so we have determined 
the maximum and minimum temperature greybodies consistent with this mid-infrared upper 
limit and the 1.1 mm and submm detections. We find a maximum temperature of 63 K, 
and a long wavelength dust emissivity index of (5 = 0.6. The maximum luminosity is then 
<3.3 L0, the mass is 6.4 x 1O~^M0, and the bolometric temperature is <70 K. The lower 
limit to the source temperature is 23 K, with a dust emissivity index of f3 = 0.8. This 
corresponds to a minimum luminosity of 0.15 L©, and a mass of 4.5 X 10-2 Mq. The ratio of 
submillimeter luminosity to total luminosity in SMA2 is between 3.4 x 10^'^ (cold limit) and 
6.6 X 10"^ (warm limit). The near-infrared flux density upper limits are between 2.5 x 10^ 
(cold limit) and 3.3 x 10^ (warm limit) times less than the S EP peak. Both o f these ratios 
are consistent with typical characteristics of a Class source ( iFroebrich II2005I ). 



Considering the 7 mm data of IChoi fc Tatematsu I (120041 ) . we note that our SED predicts 
1.9 mJy from SMA2 at that wavelength, whereas the measured value is 2.4 ± 0.7 mJy. We 
also note that the positions agree to within 0''9, with most of the offset in declination. 
Considering the low-elevation of the source when observed at the VLA, the offset could be 
due to instrumental uncertainty. These facts leads us to favor the interpretation that much 
of the emission from CT2 is dust emission from SMA2. A lesser portion of the 7 mm flux 
density (~ 0.5 — 1 mJy) could be due to free-free emission from an ionize d jet as it would still 



be consistent with a relatively fiat centimeter-wave spectral index (e.g. Reynolds ( 19861 )) 



4.2. Non-detections of other 7mm sources 



In addition to CT2, IChoi fc Tatematsu I (120041 ) detected four more point sources in the 
IRS 7 field at 7 mm with the VLA which we do not detect with the SMA. CT4 has a negative 
spectral index, a ~ —0.8 ± 0.1, i.e. the emissioii appears non-thermal, possibly gyrosyn- 
chroton emission, although iFeigelson et al.l ( 19981) did not detect any polar ized emission. 
However, inspection of the data by iBrownl (119871 ) and IFeigelson et al.l (119981 ) suggest that 
the centimeter VLA emission is variable; it could therefore still be free-free emission if the 
source was in a low state at the time of the 7 mm observations. Due to their non-detection 
at cm wavelengths, CT 1, CT 3, and CT 5 were hypot hesized to have steep spectra (a > 2.3). 
With no IR counterparts, IChoi fc Tatematsu I (120041 ) suggest that the emission originates in 
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compact dust disks around young low-mass stars. CT 1 could be associated with a faint 450 
yum peak seen to the northeast of SMA2 (Figure [3]), but we definitely see no 450 fim excess 
from CT5. Furthermore, the SMA should have detected these objects easily if they have 
temperatures typical of dust disks (20-30K). If the 7mm emission from these objects is from 
dust, the only way their SEDs can be consistent with the SMA non-detections and the VLA 
3.6cm non-detections is if they are very cold (T ~ 6K) and remain optically thick (r > 1) 
at wavelengths shortward of 1 cm. In this scenario, they could represent compact accretion 
cores just pr ior to the formation of a stellar embyro, such as the lOK "Class -1" objects 
predicted by IBoss &: Yorkd (Il995l ). In any case, the reality and nature of these potentially 
interesting objects warrant further observations with ATCA or eVLA. 



4.3. Are SMA 1 and SMA 2 Class Protostars? 



Both iNutter et al.l (120051 ) and ISchoier et al.l (120061 ) postulate that either one or both of 
the young stell a r obj ects w hich we associa t e with SMA 1 and SMA 2 are Class protostars. 
Forbrich et al.l (120061 ) and iForbrich et al.l (120071 ) suggest that SMA 1 is a candidate Class 
source, but do not detect SMA 2 in their VLA and Chandra observations. Our data 
suggest that SMA 2 is a Class source, while SMA 1 straddles Class and Class I. It is 
clear that both SMA 1 and SMA 2 are submillimeter bright objects with centimeter radio 
counterparts, hard X-ray counterparts, and no detected near-IR counterparts. We see no 
evidence for binary association. Neither object is detected at near-infrared wavelengths, but 
SMA 2 is detected in all four IRAC bands. At 24 /im we detect SMA 1 and can only place 
an upper limit on the flux of SMA 2, but the proximity to two extremely bright sources 
(R CrA and IRS 7) have likely masked any 24 /im flux that may originate from SMA 2. 
With a bolometric temperature upper limit of 70 K, no IRAC or near-IR detections, and the 
presence of a bipolar outflow, SMA 2 is very likely a Class source. SMA 1 has a bolometric 
temperature of 83 K and shows evidence for a higher temperature (350 K) component in 
its SED. This suggests a more evolved evolutionary state, consistent with a Class 0/Class I 
transition object. Its significantly higher luminosity (I7L0) also suggests that SMA 1 is a 
more ma ssive protostar. ^y hen placed on a bolometric luminosity vs. temperature (BLT) 
diagram (IMyers et al.lll998l ). SMA 2 is consistent with the evoluti o nary track of a <O.3M0 
envelope when compared to the model tracks of lYoung fc Evans I (120051 ). SMAl, with it's 
slightly higher bolometric temperature, but far higher bolometric luminosity, is consistent 
with a ~1 - 5 Mq envelope mass. Caution must be used with the use of a BLT diagram to 
derive masses, since Tboi is not a good estimator of age, resulting in the rather wide range of 
acceptable masses for the two sources. 
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4.4. A young embedded low-mass cluster? 



Schoier et al.l (|2006[) believe that SMA 1 and SMA 2 form a protobinary system similar to 



IRAS 16293 (IWalker et al.lll986l ). We do not find evidence that supports this conclusion. The 
two sources seem to be at different evolutionary states, with very different source properties. 
Our outflow data shows evidence for one large scale outflow associated with SMA 2, and at 
least one smaller outflow (associated with HHIOO-IR). A small outflow with an axis roughly 
orthgonal to the SMA 2 outflow could be driven by SMA 1. In addition, the "shell" of HCO+ 
emission seems to be spatially centered on SMA 2, with little association with SMA 1. In 
the situation of a coeval protobinary, any shells of material would likely be surrounding the 
entire protobinary system, not a single member. 

We postulate that the HCO^ emission could be the result of the interaction of many 
molecular outflows, or could trace a shell of material around SMA 2, swept up by the outflow 
or by a lower velocity wind from a single object. Would such a shell be gravitationally 
bound to the central object, or would it be expanding? Using the mass estimate from the 
previous section, we can compare the kinetic energy of the possible shell to its gravitational 
potential energy. Here, we have maximized the mass estimate of the shell, and use the 
best estimates for the protostellar envelope mass and protostar mass derived from the BLT 
diagram. We calculate that the kinetic energy is more than a factor of 20 higher than the 
gravitational potential energy. This suggests that the material is not gravitationally bound 
to the central object. Further understanding of the kinematics and structure of the dense 
molecular gas in this region will require interferometric observations combined with zero 
spacing data from properly- sized single dish telescopes, a technique that is envisioned a s 



a standard capability of the Atacama Large Millimeter Array (ALMA) (jWootten 1 120071 ). 
Similarly, mosaiced interferometric CO observations with zero spacing, or data from new, 
large single dish telescopes such as the Large Millimeter Telescope (LMT) should be able to 
more accurately trace the source (s) of the bipolar outflows. 



5. Summary 



We have observed the R CrA molecular cloud in the region surrounding the near-lR 
source IRS 7 with the Submillimeter Array in 1.1 mm continuum and the HCO"*" J = 3 ^ 2 
line. We present complementary SCUBA/ JCMT 850 /im and 450 /xm continuum images, 
Spitzer archival IRAC and MIPS 24 /im images and DesertStar/HHT ^^CO J = 3 ^ 2 
maps of the molecular outflo w. We detect two compact continuum sources, SMA 1 and 
SMA 2 previously identified bv iNutter etaD koO^ as SMMIB and SMMIC, which are not 
associated with any known near-lR counterparts but are detected at centimeter and hard 



-17- 



X-ray wavelengths. The brighter, spatially-extended source SMMIA is not detected in our 
interferometric observations. The SED of SMA 1 derived from the SMA, SCUBA and Spitzer 
data is best fit by a two component greybody with temperatures of 55 K and 344 K, and 
a long wavelength dust grain emissivity index f3 = 1.3. The bolometric temperature of the 
source is 83 K, with a bolometric luminosity of 17 ± 6 L0. While SMA 2 is detected at 1.1 
mm, 870 /xm and 450 /^.m, we do not detect it at 24 /xm or in any IRAC bands. From the 
flux density upper limit at 24 /xm, we compute the maximum temperature greybody model 
consistent with the data to be 63 K with /3 = 0.6. The bolometric temperature is <70 
K, and the bolometric luminosity is between 0.15 and 3.3 Lq. These results suggest that 
SMA 2 is a low mass (< 0.3 Mq) Class source, while SMA 1 is a transitional Class 0/Class I 
object of higher envelope mass (~l-5 Mq). We see evidence for a single large scale molecular 
outflow driven by SMA 2, with evidence for at least one more smaller outflow associated 
with HH-IOOIR (IRSl). The HCO+ J = 3 ^ 2 emission detected with the SMA is very 
clumpy in nature, with a complex velocity field consistent with a spherical shell of material 
blown by a wind from one or more of the members of the young cluster. This emission could 
also be explained by the interaction of outflows driven by the many YSOs in the region. 
Future interferometric observations with zero spacing and short baseline data will assist in 
determining the distribution of gas around these two protostars. 

C.E.G. is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship 
under award AST-0602290. We wish to thank Dr. Remy Indebetouw of the University of 
Virginia with his assistance in the PSF fitting photometry of the saturated MIPS 24 //m 
obervations. 

REFERENCES 

Andrews, S. M., & WiUiams, J. P. 2005, ApJ, 631, 1134A 
Boss, A. P., & Yorke, H. W. 1995, ApJ, 439, L55 
Brown, A. ApJ, 322, L31 

Caratti o Garatti, A., Giannini, T., Nisini, B., & Lorenzetti, D., 2006, A&A, 449, 1077 
Chen, W. P., & Graham, J. A. 1993, ApJ, 409, 319 

Chini, R., Kampgen, K., Reipurth, B., Albrecht, M., Kreysa, E., Lemke, R., Nielbock, M., 
Reichertz, L.A., Sievers, A., & Zylka, R. 2003, A&A, 409, 235 



-18- 



Choi, M. & Tatematsu, K. 2004, ApJ, 600, L55 

Cohen, M., Harvey, P. M., Wilk ing, B. A., & Schwartz, R. D., 1984, ApJ, 281, 250 

Davis, C. J., Smith, M. D., Eisloffel, J., & Davies, J. K., 1999, MNRAS, 308, 539 

Feigelson, E.D., Carkner, L., & Wilking, B.A. 1998, ApJ, 494, 215 

Forbrich, J., Preibisch, T., & Menten, K. M. 2006, A&A, 446, 155 

Forbrich, J., et al. 2007, A&A, 464, 1003 

Froebrich, D. 2005, ApJS, 156, 169 

Graham, J. A. 1993, PASP, 105, 561 

Groppi, C.E., Kulesa, C, Walker, C., & Martin, C.L. 2004, ApJ, 612, 946 

Groppi, C.E., Walker, C.K., Kulesa, C., Golish, D., Hedden, A., Narayanan, G., Lichten- 
berger, A.W., Kooi, J.W., Graf, U.U., Heyminck, S., 2004 Proc. SPIE, 5498, 290 

Hamaguchi, K., Corcoran, M.F., Petre, R., White, N.E., Stelzer, B., Nedachi, K., Kobayashi, 
N., & Tokunaga, A.T. 2005, ApJ, 623, 291 

Harju, J., Haikala, L.K., Mattila, K., Mauersberger, R., Booth, R.S., & North, H.L. 1993, 
A&A, 278, 569 

Hartigan, P., & Graham, J. A. 1997, AJ, 93, 913 

Henning, T., Launhardt, R., Steinacker, J., & Thamm, E. 1994, A&A, 291, 546 
Hildebrand, R.H. 1983, QJRAS, 24, 267 

Hogerheijde, M.R. SFChem 2002 Eds. Charles L. Curry & Michel Fich. 2003, NRG Press, 
Ottawa, Canada, p. 44 

Holland, W.S., Robson, E.I., Gear, W.K., Cunningham, C, Lightfoot, J.F., Jcnncss, T., 
Ivison, R.J., Stevens, J.A., Ade, P.A.R., Griffin, M.J., Duncan, W.D., Murphy, A., & 
Naylor, D.A. 1999, MNRAS, 303, 659 

Jenness, T., & Lightfoot, J.F. 1999, Starlink User Note 216.6, Rutherford Appleton Labo- 
ratory, Particle Physics & Astronomy Research Council 

Knude, J., & H0g, E., 1998, A&A, 338, 897 



-19- 



Koyama, K., Hamaguchi, K., Ueno, S., Kobayashi, N., & Feigelson, E. D. 1996, PASJ, 48, 
L87 

Levreault, R.M. 1988, ApJS, 67, 283 
Loren, R.B. 1979, ApJ, 227, 832 

Megeath, S.T., Allen, L.E., Gutermuth, R.A., Pipher, J.L., Myers, P.C., Calvet, N., Hart- 
mann, L., MuzeroUe, J., Fazio, G.G. 2004, ApJS, 154, 367 

Myers, P.C., & Ladd, E.F., 1993, ApJ, 413, L47 

Myers, P.C., Adams, F.C., Chen, H., & Schaff, E. 1998, ApJ, 492, 703 
Nutter, D.J., Ward-Thompson, D., & Andre, P. 2005, MNRAS, 357, 975 
Reynolds, S. P. 1986, ApJ, 304, 713 

SandeU, G., Jessop, N., & Jenness, T. 2001, "The SCUBA map reduction Cookbook", Star- 
hnk Cookbook 11.2 , Rutherford Appleton Laboratory, Particle Physics & Astronomy 
Research Council 

Sault, R.J., Teuben, P.J., & Wright, M.C.H., 1995, in ASP Conf. Ser. 77: Astronomical 
Data Analysis Software and Systems IV, Eds. R.A. Shaw, H.E. Payne, & J.J.E. Hayes 
(Astronomical Society of the Pacific: San Francisco), p. 433 

Schoier, F.L., J0rgensen, J.K., Pontoppidan, K.M., & Lundgren, A. A. 2006, A&A, 454, L67 

Snell, R.L., Scoville, N.Z., Sanders, D.B., & Erickson, N.R. 1984, ApJ, 284, 176 

Strom, S.E., Grasdalen, G.L. & Strom K.M., 1974, ApJ, 191, 111 

Taylor, K.N.R., & Storey, J.W.V. 1984, MNRAS, 209, 5 

van den Ancker, M., 1999, Ph.D. Thesis, University of Amsterdam 

Walker, C.K., Lada, C.J., Young, E.T., Maloney, P.R., & Wilking, B.A. 1986, ApJ, 309, L47 
Wang, H., Mundt, R., Henning, T., Apai, D. 2004, ApJ, 617, 1191 

Wilking, B. A., McCaughrean, M. J., Burton, M. G., Giblin, T., Rayner, J. T., & Zinnecker, 
H. 1997, AJ, 114, 2029 

Wilking, B. A., Taylor, K. N. R., & Storey, J. W. V. 1986, AJ, 92, 103 



- 20 - 

Wootten, A. 2007, ArXiv Astrophysics e-prints, arXiv:astro-ph/070266^ 
Young, C.H. & Evans, N.J. 2005, ApJ, 627, 293 

Zhao,J-H. 2006, "Miriad User's Guide to the Reduction of SMA Data" 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



- 21 - 



Table 1. SMA 271 GHz continuum sources 



Source 


q;(2000.0) 

I'/i m sj 


5(2000.0) 
[° ' "] 


ea 

[" 


X Ob P.A. 

X "] [°] 


da X dh 

[AU X AU] 


S(l.l mm) ^ 
[mJy] 


SMAl 
SMA 2 


19 01 56.416 
19 01 55.289 


-36 57 27.85 
-36 57 17.03 


1.1 


X 0.5 52 
<1 


220 X 140 
<170 


300 ± 17 
280 ± 17 



^The quoted uncertainty does not include the 20% uncertainty in the overall flux 
scale. 



Table 2. Positions, deconvolved sizes and integrated flux den sities of SCUBA so urces in 
the R Cr A cloud core following naming convention in ( INutter et al.l 120051 ) 



Sub-mm 


a(2000.0) 


.5(2000.0) 


da 


X 




P.A. 


S{850 fim) 


S(450 tim) 


Other Designations 


source 


j^/i m s j 




[" 


X 


"] 


[° ] 


[Jy] 


[Jy] 




SMM 1 A^' 


19 01 55.23 


-36 57 46.0 


43 


X 


27 


-80 


14.5 ± 2.2 


151 ± 30 


vdA 2&4 


SMM IAS 


19 05 54.97 


-36 58 29.0 


14 


X 


13 




1.2 ± 0.2 


6.9 ± 1.7 




SMM IB'' 


19 01 56.35 


-36 57 28.1 


3.1 


X 


1.3 




1.5 ± 0.2 


5.2 ± 0.5 


vdA5, SMAl, 14, BlOE, 


SMMIC 


19 01 55.29 


-36 57 17.5 


3 


X 


<1 


-9 


1.1 ± 0.2 


2.4 ± 0.7 


vdA3, SMA2, B9, FCW6, CT2&3 


SMM 2 


19 01 58.57 


-36 57 09.4 


7 


X 


6 




1.2 ± 0.2 


5.8 ± 1.5 


vdA6, 13, WMB55 


SMM 3^= 


19 01 50.81 


-36 58 11.3 


9 


X 


7 




1.5 ± 0.2 


>5.4 ± 1.0 


vdAl, IRSi, HHIOO-IR, I7,B8, FCW4, K6 


SMM 4 


19 01 47.97 


-36 57 17.6 


40 


X 


29 


-4 


4.8 ± 0.7 




IRS 5, 110, 19 



''Resolved into two cores at high angular resolution, see e.g. Fig. 2 and van den Ancker (1999) 

''Appears extended at 450 /im consistent with SMA results, but the errors from the Gaussian fits are large 

■^At the edge of the 450 fim map; Coordinates & size taken from 850 fim, flux density at 450 fim is underestimated. 



Note. — The prefixes used for the source n umbers i n column 'Other Des ignati ons' come from: vdA llvan den Ancker Ill999h. SMA CT able 
1), I (IRAQ sources , Ta ble 3). B llBrownlll987l). FCW llFeieelson et al.lll99d) . CT llChoi fc Tatematsu |[2004h . WMB dWilking et al]|l997h . K 
llKovama etal.lll996l) . X llHamaguchi et ahlboosl) 
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Table 3. Positions and flux densities of IRAC 8 fim. sources in the R CrA cloud core. 

Photometric uncertainties are ~10%. 



nr 


Optical /near-IR 


a(2000.0) 


5(2000.0) 


S(3.6 fj.m) 


S(4.5 Atm) 


S(5.8 Aim) 


S(8 Aim) 


comment 




counterpart 


[h m s] 


[° ' "] 


[mjy] 


[mJy] 


[mJy] 


[mJy] 




1 


T CrA 


19 01 58.77 


-36 57 50.0 




2636.1 


2772.3 


3113.8 


saturated at 3.6 /im 


2 


IRS 4 


19 01 58.64 


-36 56 17.1 


77.9 


55.3 


45.9 


22.2 




3 


WMB 55 


19 01 58.54 


-36 57 08.7 


26.7 


71.0 


143.2 


189.2 




4 




19 01 56.39 


-36 57 28.2 


41.9 


251.1 


687.3 


1034.2 




5^ 


IRS 7 


19 01 55.32 


-36 57 22.1 


58.2 


245.5 


666.4 


1099.3 




6 


R CrA^ 


19 01 53.65 
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^Photometry more uncertain, because the source emission is affected by one of the diffraction spikes from R CrA 
''Position from 3.6 urn, one saturated pixel; optical position (Simbad): RA = 19'' 01"" 53f 65, Dec = -36° 57' 07'.'6 
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Fig. 1. — SMA continuum interferometric image at 271 GHz. Two sources are detected 
of approximately equal flux density (0.3 Jy and 0.28 Jy). Crosses mark the positions of 
centimeter detections from iFeigelson et al.l ( 1l998l ).Tl ie open circle mark s the position of 
IRS 7, while the boxes mark the 7 mm sources from IChoi fc Tatematsu I ( 120041 ). The flUed 
circle marks the location of the IRAC IR source associated with SMA 1. Contours are spaced 
at 4cr, with the lowest contour at 4cr (36 mJy). 
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Fig. 2.— 850 fj,m (Left) and 450 //m (Right) SCUBA maps of the R CrA cloud core. Star 
symbols mark the positions of IRAC 8 /xm sources and the crosses show the two SMA sources. 
These two sources as well as IRS 7 and the HAEBE star R CrA are labeled on the 450 fim 
image. Only the 450 fim data clearly resolves SMAl and SMA 2. The minimum contour for 
the 850 fira image is 0.15 Jy, with contours separated by 0.2 Jy. For the 450 /^m image, the 
first two contours are 1.4 and 2.5 Jy, with other contours spaced by 0.9 Jy. 
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Fig. 3. — Blow-up of the SCUBA 450 /xm continuum map (grayscale and contours) showing 
SMAl, IRS 7 and SMA2. Crosses mar k the SMA positions, stars the IRAC positions, filled 
squares th e centimeter VLA position s (IFeigelson et al.lll998l ). filled triangles the 7 mm VLA 
positions (IChoi &: Tatematsu II2004I ). The first two contours are 1.4 and 2.5 Jy, with other 
contours spaced by 0.9 Jy. 
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Right Ascension (2000.0) 

Fig. 4. — IRAC 4.5 /im image of the R CrA region corresponding to figure [2l This image 
is a composite of both long and short exposure time data to improve dynamic range. A 
logarithmic stretch has been applied to emphasize fainter features. Sources are plotted 
identically to figure [21 



- 28 - 



150 
100 
50 

-50 
-100 
-150 



1 1 1 1 1 1 1 1 \j 1 1 1 1 1 1 
: +HH99B \ 


1 1 1 1 1 1 1 1 1 1 1 1 1 i_ 


+ HH99A \ 

\ 

- \ 


- 
- 






-_ """^"^"^"^^i^^^^p 


^"Tly / HH104-C_t____: 
\ HH104-D- 






( . ^ 
\ / 




- Beam 




: • 


im97 J 


-1 < 1 1 1 1 1 1 1 < 1 < 1 1 1 1 1 1 < 


, 1 , , , , 1 , , , \ , ,- 



150 100 50 -50 -100 -150 
RA Offset [Arcseconds] 

Fig. 5. — ^^CO J = 3 ^ 2 outflow map of R CrA. The map covers an area of 5' x 5'. The 
blue-shifted emission is integrated over the velocity interval —10 to 3 km s~^, while the red 
interval is 10 to 21 km s~^. Contours are separated by two sigma, with the lowest contour 
at two sigma, calculated from the noise in the spectra off the line integrated over the same 
velocity width. Optically detected Herbig-Haro objects are marked with plus symbols, while 
IR and sub-mm sources are marked with squares. Lines connect HH objects we believe are 
related to molecular outflow lobes. 
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Fig. 6. — HCO^ J = 3 ^ 2 emission overlaid on the 271 GHz continuum contours (solid 
black). Red solid lines are integrated from 5.7 to 9.3 km s^^, and the dashed blue lines are 
integrated from 3.8 to 5.7 km s^^. Contours are separated by la, with the lowest contour 
at Icr. The HCO"'" clumps seem to be distributed roughly spherically around SMA2, and 
show no direct association with the continuum sources nor with the large scale CO outflows 
shown in Figured 
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Fig. 7. — Spectral energy distribution of the source SMA 1. 24 fim MIPS and IRAC data 
points are derived from Spitzer archival data, while the 1.1 mm, 850 /im and 450 fim data are 
from this work. The 850 /im measurement is likely influenced by emission from SMA 2 which 
is on ly barely resolved fr o m SM A 1 in the single-dish data. The 4a upper limit at 7 mm is 
from lChoi fc Tatematsu I (120041 ). The data are fit with a two temperature diluted blackbody 
plus a power-law free-free component. The upper limits at 60 and lOO/im are from the PSC 
detection of IRAS 18585-3701. They are drawn as upper limits because SMAl and SMA 2 
are unresolved in these data. Likewise, the detec tions at 52 and 100 /xm are from the Kuiper 
Airborne Observatory with a 38" to 47" beam ( ICohen et al.lll984l ). The spectral index of 
the free-free component plotted i s F,y ^ as determined by the two me a surem ents at 



centimeter wavelengths (6 cm by iBrownl (119871 ) and 3.6 cm by lForbrich et al.l (120071 )) 
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Fig. 8. — Spectral energy distribution of the source SMA2. The 850 /xm measurement is 
hkely influenced by emission from SMA 1 which is only ba rely resolved from SMA2 in the 
single-dish data. The detection at 7 mm is from source 2 of lChoi &: Tatematsu I (120041). The 
flux density at 3.6 cm is the average of ten measurements from iForbrich et al.l (120061 ) and 
Feigelson et al.l (119981 ). The data are fit with a power-law free-free component [Fu cx i>~^'^) 
plus a single temperature greybody, using the 24 /im upper limit to find the maximum 
temperature diluted blackbody consistent with the data. The spectral index of the free-free 



component plotted is F^, 



-0.18 



, the shallowest value consistent with all the data. 



